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Abstract. We isolated and characterized i 6 new di- and tetranucleotide microsatellite markers for the critically endangered 
Asian hox turtle genus Cuora, focusing on the “Cuora trifasciata” species complex. The new markers were then used to 
analyse genetic variability and divergence amongst five described species within this complex, namely C. aurocapitata (n = 
i 8 ), C. cyclornata (n = 31 ), C. pani (n = 6 ), C. trifasciata (n = 58 ), and C. zhoui (n = 7 ). Our results support the view that 
all five species represent valid taxa. Within two species (C. trifasciata and C. cyclornata), two distinct morphotypes were 
corroborated hy microsatellite divergence. For three individuals, morphologically identified as being of hybrid origin, the 
hybrid status was confirmed by our genetic analysis. Our results confirm the controversial species {Cuora aurocapitata, 
C. cyclornata) and subspecies/morphotypes (C. cyclornata meieri, C. trifasciata cf trifasciata) to be genetically distinct, 
which has critical implications for conservation strategies. 
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Introduction 

Asian box turtles of the genus Cuora are small to medium¬ 
sized species (11-35 cm straight carapace length) that in¬ 
habit rainforests in South and East Asia (Fig. 1 shows the 
morphology and distribution of the species investigated 
in this paper). Habitat destruction and collection for con¬ 
sumption or as pets are the primary threats to these semi- 
aquatic freshwater turtles (Turtle Conservation Coalition 
2011) and all but one species (C. amboinensis) are now 
listed as ‘critically endangered’ in the lUCN Red List of 
Threatened Species (lUCN 2012). Two species, C. mccordi 
and C. zhoui, might already have become extirpated in the 
wild (Turtle Conservation Coalition 2011). 

Up to 13 species have been described, but taxonomic un¬ 
certainty persists, hampering conservation efforts. Studies 
on mitochondrial DNA (mtDNA) provided a first taxo¬ 
nomic framework (e.g., Stuart & Parham 2004, Spinks 


et al. 2004, Parham et al. 2004), but later work using a few 
nuclear loci suggested that some putative species could be 
hybrids (Stuart & Parham 2007). High-resolution nucle¬ 
ar markers (like microsatellites) were not previously avail¬ 
able for the genus, however, even though these are essen¬ 
tial tools for species delimitation and the identification of 
hybrids, both of which are paramount for conservation 
projects and legislation. It is also important to develop ana¬ 
lytical tools by which confiscated turtles or animals in cap¬ 
tive breeding programs can be traced back to their popula¬ 
tion of origin. 

Taxonomic uncertainty has plagued conservation ef¬ 
forts, especially in the C. trifasciata species complex. Mor¬ 
phological and mtDNA sequence divergence suggest the 
existence of two species within what was once consid¬ 
ered Cuora trifasciata (Blanck et al. 2006): C. trifasciata 
sensu stricto (‘clade A) and C. cyclornata (‘clade C’), but 
the validity of C. cyclornata has been challenged based on 
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mtDNA and nuclear single copy loci sequences (Spinks 
& Shaffer 2007, Spinks et al. 2009, Spinks et al. 2012). 
Within C. cyclornata, Blanck et al. (2006) detected two 
divergent lineages, which they described as the subspe¬ 
cies C. c. cyclornata and C. c. meieri. Furthermore, several 
specimens of unknown origin formed yet another geneti¬ 
cally divergent C. trifasciata clade (clade B’), but clades A 
and B were indistinguishable based on the low-resolution 
nuclear markers used thus far (Spinks & Shaffer 2007). 
Even the most recent study (Spinks et al. 2012), which in¬ 
cluded sequences of 15 nuclear single-copy loci, could not 
resolve these issues. 

Yet another level of taxonomic uncertainty concerns 
C. aurocapitata and C. pani, which - although distinguish¬ 
able by morphology and some nuclear loci (Spinks et al. 
2012) - show extensive mt-haplotype sharing (Parham et 
al. 2004, Spinks et al. 2004). This suggests that they inter¬ 
breed in secondary contact zones, e.g., in western Anhui, 
southern Henan or eastern Hubei Province, from where a 
phenotypically intermediate variety is known (Blanck & 
Tang 2005, T. Blanck pers. obs.). 

Microsatellite analysis is an elegant tool to assess the 
degree of reproductive isolation, delimit taxonomic units, 
and detect potential hybridisation between taxa. The aims 
of this study were (1) to develop a new informative set of 
highly polymorphic, biparentally inherited molecular 
markers (i.e., microsatellites) for Cuora spp., and (2) to 
evaluate the degree of reproductive isolation among the 
distinguishable morphotypes of the C. trifasciata species 
complex. Here, we perform a genetic evaluation of the pre¬ 
diction of the Biological Species Concept (BSC) among 
the a priori defined morphotypes, i.e., should diagnosable 


morphotypes occur in sym- or parapatry, then significant 
genetic divergence among these morphotypes indicates re¬ 
productive isolation. We are aware of the fact that genet¬ 
ic distinctiveness between allopatric populations does not 
necessarily preclude that these populations could in prin¬ 
ciple interbreed (Milinkovitch et al. 2001 ). 

Material and methods 

We obtained blood and tissue samples from 120 Cuora 
specimens, morphologically assigned to 5 focal species, i.e., 
C. aurocapitata (n = 18), C. cyclornata (n = 31), C. pani (n = 
6), C. trifasciata (n = 58), and C. zhoui (n = 7; Fig. 1). Our 
samples included the holotypes of C. trifasciata and C. cycl¬ 
ornata (see Supplementary material). 

Genomic DNA was extracted using the DNeasy Tis¬ 
sue Kit (Qiagen). Four microsateUite-enriched genomic 
DNA libraries were constructed from a total DNA extract 
of a C. trifasciata sample following Paulus & Tiedemann 
(2003): genomic DNA was simultaneously restricted with 
Nhel, Haell, Rsal, and Haelll. After treatment with mung 
bean exonuclease and calf intestinal phosphatase, genomic 
DNA fragments were blunt-end legated to SNX-linkers and 
PCR-amplified with linker primers. Fragments were hybrid¬ 
ised with 5'-biotin-labelled microsatellite probes ([GA]^, 
[GT]^^, [AGTG]^, [AGAT]^), conjugated with streptavidin- 
coated magnetic beads (Dynal), and extracted with a mag¬ 
netic device. PGR products were cloned using the TOPO 
TA Gloning Kit for Sequencing (Invitrogen). PGR products 
of recombinants were blotted onto a nylon membrane and 
again hybridised with the microsatellite probes. Positive 





Cuora aurocapitata Clade A 
Cuora aurocapitata Clade B 
Cuora pani 

Cuora trifasciata trifasciata (A) 
Cuora trifasciata cf. trifasciata (B) 

Cuora cyclornata meieri 

Cuora cyclornata cf. cyclornata 
Cuora cyclornata cyclornata 

Cuora zhoui 


^ Out-of-Range record ^ Terra typica 0 UteratureA/illage/Market record 
A Wild/Scientific record ? suspected range outline 




Figure 1. Colouration (left to right: carapace dorsal, carapace ventral, head lateral), distribution, and sample origins of the analysed 
Cuora taxa. 
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clones were detected using the Phototope-Star chemolu¬ 
minescent detection kit (New England BioLabs), purified 
with Exonuclease I and Antarctic Phosphatase (New Eng¬ 
land BioLabs), sequenced with the BigDye Terminator Cy¬ 
cle Sequencing Kit version 3.1 (Applied Biosystems), and 
analysed on an ABI Prism 3130 multicapillary automatic 
sequencer (Applied Biosystems). Primers were designed in 
microsatellite flanking regions, and 16 primer pairs showed 
reliable amplifications in all 5 focal species. 

The variability of microsatellite loci was analysed by 
PCR using 5'-fluorescence-labelled forward primers. 
25 ql-reaction volumes (buffer solution: 10 mM Tris-HCl, 
pH 9.0, 50 mM KCl, 1.5 mM MgCl^, 0.1% Triton Xioo, 
0.2 mg/ml BSA) were prepared as follows: 2.5 ql dNTP 
mrxture (2.0 mM each), 0.5 pi of each primer (10 pM), 1 pi 
DNA template (-20 ng), and 0.1 pi Taq polymerase (5 U/ 
pi, MP Biomedicals). Cycling parameters included initial 
denaturation at 94°C for 5 min, 5 cycles at 94°C for 30 s, an 
elevated locus-specific annealing temperature T-i-3°C for 
1 min, 72°C for 45 s; 35 cycles at 94°C for 30 s, T (Tab. 1) 
for 1 min, 72°C for 45 s, and a final extension at 72°C for 
10 min. Fragment sizes were determined on an ABI Prism 
3130 multicapillary automatic sequencer, using the Gene- 
Mapper V 4.0 software and an internal size standard (LIZ 
500, Applied Biosystems, 250 bp peak omitted). 

We conducted an Analysis of Molecular Variance 
(AMOVA) to test for molecular divergence among 
Cuora species and calculated pairwise values using 
ARLEQUIN (Excoffier & Lischer 2010); significance 
was evaluated based on 1,000 permutations. Locus-wise 
Fj^ values were assessed with GENALEX 6 (Peakall & 
Smouse 2006). ARLEQUIN was also used to calculate 
observed (H^^) and expected (H^) heterozygosity and to 
test for deviation from the Hardy-Weinberg equilibrium 
(HWE) for each locus in each population using Fisher’s 
exact test and the Markov Chain method (1,000 demori- 
sation steps, 100 batches, with 10,000 iterations per batch 
set). Alpha-levels of HWE and linkage disequilibrium tests 
were Bonferroni-corrected for multiple comparisons. 

A Bayesian clustering was performed using the software 
STRUCTURE v 2.3.2 (Pritchard et al. 2000). Genet¬ 
ic subdivision was evaluated by estimating the likelihood 
and sample composition of independent runs of subgroups 
(k = 2-9), assuming an admixture model, with 100,000 it¬ 
erations discarded as burn-in, and a collection period of 
900,000 iterations. To check for convergence of the Markov 
Chain parameters, 15 replicate runs for each value of k were 
performed. To estimate relationships among individuals of 
different geographic origin, we constructed a Neighbour 
Joining (NJ) tree based on individual Nei’s minimum ge¬ 
netic distances and 100 bootstraps over loci (D ; Nei 1987) 
using POPULATIONS v 1.2.28 (Langella 2001). 

Results and discussion 

Our 16 new microsatellites represent the first loci specifical¬ 
ly developed for Cuora spp. All loci amplified successfully 


throughout morphologically distinguishable taxa. Locus- 
wise allele numbers ranged from one to 22 alleles per taxon 
(Tab. 1); no locus was monomorphic in all species though, 
and no locus was consistently heterozygote-deficient across 
species. Genetic variability was high in C. trifasciata (mean 
Hg = 0.735), C. cyclornata (0.698) and C. aurocapitata 
(0.621), but low in C. zhoui (0.399) and C. pani (0.410), al¬ 
luding to differences in the remaining effective population 
size (estimates range from ‘possibly extinct in the wild’ for 
C. zhoui, fewer than 100 specimens for C. aurocapitata, to 
fewer than 1,000 specimens for C. cyclornata, C. trifasciata 
and C. pani; Turtle Gonservation Goalition 2011). The low 
genetic variability in C. pani can be explained by the ma¬ 
jority of samples of all other species stemming from wild 
specimens, while C. pani samples mostly came from cap¬ 
tive-bred animals where there was a possibility that some 
individuals could be related. Note that estimates of genetic 
diversity in C. zhoui and C. pani are associated with a high¬ 
er statistical error due to the small sample size. In the case 
of G. trifasciata, consistent heterozygote deficiency (at 14 of 
18 loci) hints at the presence of a hidden genetic substruc¬ 
ture, leading to a Wahlund effect (see below). 

The AMOVA indicated substantial genetic differentia¬ 
tion across taxa, as 27.8% of the variance was apportioned 
to between-species differences (i.e., global = 0.278). 
This divergence was also reflected by locus-wise Fg.j. val¬ 
ues between 0.323 and 0.732 - remarkably high values for 
microsatellites. Pairwise F^^ values showed highly statisti¬ 
cally significant differences amongst all species, ranging 
from 0.112 (C. trifasciata vs. C. cyclornata) to 0.628 (C. pani 
vs. C. zhoui; Tab. 2). STRUGTURE retrieved k = 5 and k = 7 
to have the highest likelihoods; these partitions are in good 
agreement with the morphological assignment (Fig. 2). For 
k = 5, five morphospecies (C. aurocapitata, C. cyclornata, 
C.pani, C. trifasciata, and C. zhoui) formed well-supported 
clades to which most specimens were assignable with high 
likelihoods of 80 to 100%, including the holotypes of C. tri¬ 
fasciata and C. cyclornata (Fig. 2). 

While our findings are generally in agreement with the 
results of previous studies of the genus Cuora (Spinks & 
Shaffer 2007, Spinks et al. 2009, 2012), they provide ad¬ 
ditional resolution for the morphospecies C. cyclornata 
(green in Fig. 2) and C. trifasciata (orange), which were 
clearly supported as distinct units by our microsatellite 
analysis, while previous studies concluded them to be ge¬ 
netically indistinguishable. This is perfectly congruent with 
their geographic isolation (Fig. 1) and validates the mor¬ 
phological species assignment by Blanck et al. (2006). 
Only a few samples represent putative hybrids; all of them 
were of captive origin (Fig. 2). 

If we evaluate the results of the STRUCTURE analysis 
for k = 7 (Fig. 2), further genetic structure within the two 
major clades becomes apparent: within C. cyclornata, two 
of the three morphologically identified subclusters match 
the subspecies C. c. meieri and C. c. cyclornata described 
by Blanck et al. (2006) and include the respective holo¬ 
types, although there is a mismatch for a single individual 
between morphological (C. c. cyclornata) and genetic (C. c. 
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T (°C) - PCR annealing temperature; - expected heterozygosity; - observed heterozygosity; P - probability of deviation from Hardy-Weinberg Equilibrium. 
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Table 2. Genetic divergence among the five Cuora species investigated. Above diagonal: Fj.j, values; below diagonal: associated P values. 


C. aurocapitata 

C. pani 

C. trifasciata 

C. cyclornata 

C. zhoui 

C. aurocapitata 

- 

0.355 

0.290 

0.357 

0.486 

C. pani 

< 0.001 

- 

0.356 

0.404 

0.628 

C. trifasciata 

< 0.001 

< 0.001 

- 

0.112 

0.367 

C. cyclornata 

< 0.001 

< 0.001 

< 0.001 

- 

0.425 

C. zhoui 

< 0.001 

0.002 

< 0.001 

< 0.001 

- 


meieri) assignments. An additional morphological group 
(C. cyclornata cf. cyclornata) was genetically not supported. 
Likewise, two subclusters can be seen within C. trifasciata, 
one being composed of Chinese mainland specimens in¬ 
cluding the holotype, the other of animals from Hainan 
Island, which can be assigned to C. trifasciata mtDNA 
‘Clade B’ (C. trifasciata cf trifasciata; Blanck et al. 2006; 
Fig. 2). This substructure within C. trifasciata explains the 


heterozygote deficiency observed when heterozygosity is 
estimated species-wide (Wahlund effect). While the sub¬ 
clusters of C. cyclornata are easily distinguished based on 
morphological characters, morphological identification is 
far less straightforward in the cases of the two C. trifasciata 
subclusters. Further studies are needed to better character¬ 
ize the potentially new subspecies uncovered in this study. 
All 7 clusters were confirmed in a Neighbour Joining (NJ) 


K = 5 


C. aurocapitata C. pani C. U. C. cyclornata C. zhoui 



C.; L, - .If c , ■ C.l. cf. trifasciata Clade B 

C. aurocapitata C. panI C.c. meieri C.c. cyclornata C. zhoui 



K = 7 

Figure 2. Proportional assignment of Cuora specimens to genetically inferred clusters (STRUCTURE results for k = 5 [upper] and k = 
7 [lower]). Vertical black bars separate morphologically assigned species/subspecies; * depicts individuals that were morphologically 
identified as hybrids and # depicts type specimens. 
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• C. aurocapitata 
C.pani 

A C. t. trifasciata Clade A 
A C. t. cf. trifasciata Clade B 
■ C. c. meieri 
Q C. c. cyclornata 
O C. zhoui 
♦♦♦ Hybrids 



Figure 3. Individually based neighbour-joining tree of NeTs minimum genetic distance for the microsatellite data. Each individual is 
represented by a symbol corresponding to its morphological species/subspecies assignment as shown in the legend. 


tree based on individual NeTs minimum genetic distances 
(Fig. 3 )- 

We detected only one presumed subspecies hybrid with¬ 
in C. cyclornata (C. c. cyclornata x C. c. meieri), and two 
within C. trifasciata, equalling only 5% of the 92 specimens 
of these two morphospecies analysed herein. We are thus 
inclined to argue that C. cyclornata is indeed a valid species 
with a least two well-distinguishable subspecies and that 
neither C. cyclornata nor C. aurocapitata is of hybrid origin 
as has been hypothesized by Spinks et al. (2012). 

In summary, our analysis of 16 new microsatellite loci 
demonstrated that the previous morphospecies/subspecies 
assignment in Cuora spp. (e.g., Blanck et al. 2006) is ful¬ 
ly justified for the taxa examined here. These markers can 
be used, e.g. to genotype confiscated animals, and should 
form the backbone of future captive breeding initiatives, 
particularly for taxa that are morphologically difficult to 
distinguish. Also, the application of this marker system will 
help to avoid the accidental inclusion of hybrid specimens 
in captive breeding programs. We highlight the validity of 


the species status for C. cyclornata as demonstrated in this 
study and wish to draw attention to the precarious status of 
the few remaining natural C. cyclornata and C. trifasciata 
populations (Blanck et al. 2006, Turtle Conservation 
Consortium 2011). Our new analytical tools hopefully will 
prove useful in future reintroduction programs by allowing 
animals to be genotypified prior to their captive breeding 
or release. 
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SI. Information about analysed specimens. 


Taxa Taxon according to 

STRUCTURE 


European Studbook 
Foundation No. 


Additional field name/information 


1 A and B 

1 

Cuora aurocapitata 

ESF-AUR-HH 

CA Herschel/Switzerland; wUd-caught 

2 

2 

Cuora aurocapitata 

ESF-AUR-1 

CA MSI; wild-caught 

3 

3 

Cuora aurocapitata 

ESF-AUR-2 

CA MS2; wild-caught 

4 

4 

Cuora aurocapitata 

ESF-AUR-3 

CA MS3; wild-caught 

5 

5 

Cuora aurocapitata 

ESF-AUR-4 

CA MS4; wild-caught 

6 

6 

Cuora aurocapitata 

ESF-AUR-5 

CA MSS; wild-caught 

7 

7 

Cuora aurocapitata 

ESF-AUR-6 

CA MS6; wild-caught 

8 

8 

Cuora aurocapitata 

ESF-AUR-7 

CA MS7; wild-caught 

9 

9 

Cuora aurocapitata 

ESF-AUR-8 

CA MS8; wild-caught 

10 

10 

Cuora aurocapitata 

ESF-AUR-9 

CA MS9; wild-caught 

11 

11 

Cuora aurocapitata 

ESF-AUR-10 

CA MS 10; wild-caught 

12 

12 

Cuora aurocapitata 

ESF-AUR-11 

CA MSll; wild-caught 

13 

13 

Cuora aurocapitata 

ESF-AUR-12 

CA MS 12; wild-caught 

14 

14 

Cuora aurocapitata 

ESF-AUR-13 

CA MS13; wild-caught 

15 

15 

Cuora aurocapitata 

ESF-AUR-15 

CA MS15; captive-bred 

16 

16 

Cuora aurocapitata 

ESF-AUR-20 

CA MS20; captive-bred 

17 

17 

Cuora aurocapitata 

ESF-AUR-24 

CA MS24; captive-bred 

18 

18 

Cuora aurocapitata 

ZT-JIl 

C.A.l, Jingxian, Anhui; wild-caught 


19 C 

20 
21 
22 

23 

24 



ZT-PAl 
ESP-PAN-1 
ESF-PAN-2 
ESF-PAN-3 
ESF-PAN-5 
ESF-PAN-9 


PANE wUd-caught 
Cp MSI; wild-caught 
Cp MS2; wild-caught 
Cp MS3; captive-bred 
Cp MSS; captive-bred 
Cp MS9; captive-bred 


25 D 

1 

C- t. trifasciata Clade A 

OUM8557 

OUM 8557 Cuora fasciata exchno 2009001, Luofoshan, 





Guangdong, China; wild-caught 

26 

2 

G t. trifasciata Clade A 

ESF-TRI-1 

Rotterdam 701253 

27 

3 

G t. trifasciata Clade A 

ESF-TRI-2 

Rotterdam 703160 male; wUd-caught 

28 

4 

C. t. trifasciata Clade A 

ESF-TRI-3 

Cuora trifasciata 703875; wUd-caught 

29 

5 

G t. trifasciata Clade A 

ESF-TRI-7 

Cuora trifasciata 704309; wild-caught 

30 

6 

G t. trifasciata Clade A 

MNHN9102 

MNHN 9102, Macao; wUd-caught, obtained on a market 

31 

7 

G t. trifasciata Clade A 

2006-16 

TB4 Mio / CT4; wUd-caught 

32 

8 

G t. trifasciata Clade A 

2006-13 

TB7 HK2 Hong Kong / CT7; wild-caught 

33 

9 

G t. trifasciata Clade A 

2006-19 

CT13 (Bill E/ 19); wUd-caught 

34 

10 

G t. trifasciata Clade A 

2006-20 

CT14 (BiU H / 20); wild-caught 

35 

11 

G t. trifasciata Clade A 

2006-10 

CT17 (BiU L / 17); wUd-caught 

36 

12 

G t. trifasciata Clade A 

2006-11 

CT19 (BiU N / 15); wild-caught 

37 

13 

G t. trifasciata Clade A 

2006-12 

CT20 (BiU 0 / 16); wild-caught 

38 

14 

G t. trifasciata Clade A 

2006-18 

TB15 NMW Canton / CT26, Guangzhou, China; 





wild-caught, obtained on a market 

39 

15 

G t. trifasciata Clade A 

ESF-TRI-12 

CT MS8 (=CM8 Poschadel); wUd-caught 

40 

16 

G t. trifasciata Clade A 

ESF-TRI-13 

CT MS13 (=CM13 Poschadel); wUd-caught 

41 

17 

G t. trifasciata Clade A 

ESF-TRI-42 

Pinky; wild-caught 

42 

18 

G t. trifasciata Clade A 

ESF-TRI-43 

Perky; wild-caught 

43 

19 

G t. trifasciata Clade A 

ESF-TRI-61 

wild-caught 

44 

20 

G t. trifasciata Clade A 

ESF-TRI-62 

wild-caught 

45 

21 

G t. trifasciata Clade A 

ESF-TRI-63 

wild-caught 

46 

22 

G t. trifasciata Clade A 

ESF-TRI-64 

wild-caught 

47 

23 

G t. tr^sciata Clade A 

ESF-TRI-65 

wild-caught 
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Taxa 


Taxon according to 
STRUCTURE 

European Studbook Additional field name/information 

Foundation No. 

48 D 

24 

C t. trifasciata Clade A 

ESF-TRI-300 

wild-caught 

49 

25 

C t. trifasciata Clade A 

ESF-TRl-57 

wUd-caught 

50 

26 

C t. trifasciata Clade A 

ESF-TRI-56 

wild-caught 

51 

27 

C t. trifasciata Clade A 

ESF-TRl-59 

wUd-caught 

52 

28 

C t. trifasciata Clade A 

ESF-TRI-58 

wild-caught 

53 

29 

C t. trifasciata Clade A 

ESF-TRl-60 

wUd-caught 

54 

30 

C. t. trifasciata Clade A 

ESF-TRI-49 

wild-caught 

55 

31 

C t. trifasciata Clade A 

ESF-TRl-344 

wUd-caught 

56 

32 

Cl t. trifasciata Clade A 

ESF-TRI-100 

wild-caught 

57 

33 

C. t. trifasciata Clade A 

ESF-TRI-112 

wild-caught 

58 

34 

C. t. trifasciata Clade A 

ESF-TRl-123 

captive-bred 

59 

35 

C t. trifasciata Clade A 

ESF-TRI-124 

captive-bred 

60 

36 

C t. trifasciata Clade A 

ESF-TRl-113 

wUd-caught 

61 

37 

C. t. trifasciata Clade A 

ESF-TRI-114 

wild-caught 

62 

38 

C t. trifasciata Clade A 

ESF-TRl-115 

wUd-caught 

63 

39 

C t. trifasciata Clade A 

ESF-TRI-116 

wild-caught 

64 

40 

C t. trifasciata Clade A 

ESF-TRl-117 

wUd-caught 

65 

41 

S. t. trifasciata Clade A 

ESF-TRl-129 

captive-bred 

66 

42 

C t. trifasciata Clade A 

ESF-TRI-134 

captive-bred 

67 

43 

C t. trifasciata Clade A 

ESF-TRl-81 

wUd-caught 

68 

44 

C. t. trifasciata Clade A 

ESF-TRI-47 

wild-caught 

69 

45 

S t. trifasciata Clade A 

ESF-TRl-219 

captive-bred 

70 


Morph, hybrid 

ESF-TRI-220 

captive-bred 

71 


Morph, hybrid 

JLl 

captive-bred 


72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

1 

2 

3 

4 

5 

6 

7 

8 
1 
2 

3 

4 

5 

6 


C. c. meieri 
C. c. meieri 
C. c. meieri 
C c. meieri 
C. c. meieri 
C. c. meieri 
C. c. meieri 
C c. meieri 
Cuora cyclornata cf. 
Cuora cyclornata cf. 
Cuora cyclornata cf 
Cuora cyclornata cf 
Cuora cyclornata cf 
Cuora cyclornata cf 



cyclornata 

cyclornata 

cyclornata 

cyclornata 

cyclornata 

cyclornata 


ESF-TRI-4 Cuora trifasciata 703564; wild-caught 

ESF-TRI-6 Cuora trifasciata 704308; wUd-caught 

MNHN9103 MNHN 9103, Macao; wUd-caught, obtained on a market 
2006-22 CT15 (Bill J / 22); wild-caught 

2006-23 CT16 (BUI K/ 21); wUd-caught 

2006-21 CT18 (Bill M / 23); wild-caught 

ESF-TRI-9 CT MSIO (=CM10 Poschadel); wUd-caught 

ESF-TRI-10 CTMSll (=CM11 Poschadel); wild-caught 

ESF-TRI-11 CT MS12 (=CM12 Poschadel); wUd-caught 

ESF-TRI-37 Chavornay; wild-caught 

ESF-TRI-48 captive-bred 

ESF-TRI-44 wUd-caught 

ESF-TRI-111 wild-caught 

2006-8 CT23 (BUI S / 8); wUd-caught 

2006-9 CT24 (Bill U / 9); wUd-caught 

ESF-CYC-20 CT MS4 (=CM4 Poschadel); wUd-caught 

ESF-CYC-21 CT MS5 (=CM5 Poschadel); wUd-caught 

ESF-CYC-22 CT MS6 (=CM6 Poschadel); wUd-caught 

ESF-CYC-24 CT MS29; wild-caught? 

MHNT T69 Tam Dao, Vietnam; wild-caught 

ESF-CYC-23 wUd-caught 

2006-5 TB12 Ting / CT12; captive-bred 

2006-6 CT21 (BUI P / 4); wUd-caught 

2006-4 CT22 (BUI R / 5); wild-caught 

VI Vu Quang, Vietnam; wUd-caught 

ESF-CYC-15 CT MSI (=CM1 Poschadel); wUd-caught 

ESF-CYC-16 CT MS2 (=CM2 Poschadel); wUd-caught 
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99 

G 

7 

Cuora cyclornata cf. cyclornata 

ESF-CYC-17 

CT MS3 (=CM3 Poschadel); wild-caught 

100 


8 

Cuora cyclornata cf cyclornata 

ESF-CYC-18 

CT MS7 (=CM7 Poschadel); wild-caught 

101 


9 

Cuora cyclornata cf. cyclornata 

ESF-CYC-19 

CT MS41; wild-caught 

102 


10 

Cuora cyclornata cf cyclornata 

ESF-CYC-55 

wild-caught 

103 


11 

Cuora cyclornata cf cyclornata 

ESF-CYC-53 

captive-bred 

104 


12 

Cuora cyclornata cf cyclornata 

ESF-CYC-54 

captive-bred 

105 


13 

Cuora cyclornata cf cyclornata 

ZT-CYC2 

TB5 Farm Schuppe / CT5; captive-bred 

106 

H 

1 

C. c. cyclornata 

2006-2 

“1”, Hue, Vietnam; wild-caught 

107 


2 

C. c. cyclornata 

ESF-CYC-50 

“3”, Vietnam; wild-caught 

108 


3 

C. c. cyclornata 

ESF-CYC-52 

“5”; wild-caught 

109 


4 

C. c. cyclornata 

ZFMK 71348 

TB2 ZFMK / CT2, Phong Nha Ke Bang, Vietnam; 
wild-caught 

110 


5 

C. c. cyclornata 

ZT-CYCl 

TBll Ting / CTlls; wild-caught 

111 


6 

C. c. cyclornata 

ESF-CYC-14 

CT MS9 (=CM9 Poschadel); wild-caught 

112 


7 

C. c. cyclornata 

ESF-CYC-66 

captive-bred 

113 


8 

C. c. cyclornata 

ESF-CYC-67 

captive-bred 

114 


9 

C. c. cyclornata 

ESF-CYC-68 

captive-bred 

115 


10 

C. c. cyclornata 

ESF-CYC-51 

captive-bred 

116 

I 

1 

Cuora zhoui 

ESF-ZHO-1 

CZ MSI; wUd-caught 

117 


2 

Cuora zhoui 

ESF-ZHO-2 

CZ MS2; wild-caught 

118 


3 

Cuora zhoui 

ESF-ZHO-3 

CZ MS3; wUd-caught 

119 


4 

Cuora zhoui 

ESF-ZHO-4 

CZ MS4; wild-caught 

120 


5 

Cuora zhoui 

ESF-ZHO-5 

CZ MS5; wUd-caught 

121 


6 

Cuora zhoui 

ESF-ZHO-6 

CZ MS6; wild-caught 

122 


7 

Cuora zhoui 

NTM-9001 

ZHOUI2X, Pingxiang, China; wUd-caught 


S2. Cuora taxa included in study. 


Number of 
specimens 


A Cuora aurocapitata Clade A 
B Cuora aurocapitata Clade B 

D C. t. trifasciata Clade A 45 

Morph, hybrid 2 

F Cuora cyclornata meieri 8 

G Cuora cyclornata cf. cyclornata 13 

H Cuora cyclornata cyclornata 10 


I Cuora zhoui 
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